for cc = 1:length(place)

    qualities_vec.(place{cc}) = randn(n_varieties,1);
    qualities_vec.(place{cc})(qualities_vec.(place{cc})<0) = -1*qualities_vec.(place{cc})(qualities_vec.(place{cc})<0);
    qualities_vec.(place{cc})=10000*qualities_vec.(place{cc})+1;
    qualities_vec_consumed.(place{cc})=qualities_vec.(place{cc});
    domestic_competitor_qual_vec.(place{cc})=qualities_vec.(place{cc});
    next_best_qual_vec.(place{cc}) = qualities_vec.(place{cc});

%Set probability vector & matrix
    pstolenbyentrant.(place{cc}) = eta.(place{cc});
    pstolenbyincumbent.(place{cc}) = lambda.(place{cc});
%     pkept.(place{cc}) = 1 - mu.(place{cc});
%     probvec.(place{cc}) = [pstolenbyentrant.(place{cc}) ...
%                            pstolenbyincumbent.(place{cc})...
%                            pkept.(place{cc})];
%     probmat.(place{cc}) = repmat(probvec.(place{cc}),n_varieties,1);

%Each variety has an owner (producer)
    varown_vec.(place{cc}).('y1') = [1:n_varieties]';
    varown_vec_dormant.(place{cc}) = nan(n_varieties,1);
    dormant_still_alive.(place{cc}) = zeros(n_varieties,1);
%Track the existence and size of each firm (number of firms is changing over time too);
    firm_list.(place{cc}) = [1:n_varieties]';
    firm_birthday.(place{cc}) = ones(size(firm_list.(place{cc})));
    firm_exit.(place{cc}) = nan(size(firm_list.(place{cc})));
    countvars.(place{cc}) = ones(size(varown_vec.(place{cc}).('y1')));
    innovationrate_vec.(place{cc}) = [];
    growthrate_vec.(place{cc}) = [];
    geometricgrowthrate_vec.(place{cc}) = [];

    %For tracking distributional movements
    num_active_firms.(place{cc})=[];
    num_exporters.(place{cc})=[];
    num_stopped_exporters.(place{cc})=[];
    num_non_exporters.(place{cc})=[];
    num_new_exporters.(place{cc})=[];
    num_all_new_exporters.(place{cc})=[];
    num_all_stopped_exporters.(place{cc})=[];
    num_nontraded.(place{cc})=[];
    num_exportED.(place{cc})=[];
    num_nontraded_lost_for.(place{cc})=[];
    num_nontraded_lost_dom.(place{cc})=[];
    num_exportED_lost_for.(place{cc})=[];
    num_produced_lost_for.(place{cc})=[];
    num_exportED_lost_dom.(place{cc})=[];
    qual_consumed_growth_rate.(place{cc})=[];
 
    labor_share.(place{cc})=[];
    export_share.(place{cc})=0.1017;
    mu_bar.(place{cc})=1.25;
    std_dev_firm_output.(place{cc})=[];
    agg_export_revenue_share.(place{cc})=[];
    std_dev_firm_TFPR.(place{cc})=[];
    var_firm_TFPR.(place{cc})=[];
    mean_employment_exporters.(place{cc})=[];
    mean_employment_non_exporters.(place{cc})=[];
    mean_firm_TFPR_exporters.(place{cc})=[];
    mean_firm_TFPR_non_exporters.(place{cc})=[];
    std_dev_firm_TFPR_exporters.(place{cc})=[];
    std_dev_firm_TFPR_non_exporters.(place{cc})=[];
    share_added_output_expanding_incumbents_5yr.(place{cc})=[];
    share_lost_output_contracting_incumbents_5yr.(place{cc})=[];
    share_employment_exiters_5yr.(place{cc})=[];
    share_employment_dom_entrants5yrs.(place{cc})=[];
    
    share_exporters_that_exit5yr.(place{cc})=[];
    share_exporters_that_stop_exporting5yr.(place{cc})=[];
    share_exporters_that_keep_exporting5yr.(place{cc})=[];
    share_non_exporters_that_exit5yr.(place{cc})=[];
    share_non_exporters_that_start_exporting5yr.(place{cc})=[];
    share_non_exporters_that_keep_not_exporting5yr.(place{cc})=[];
    share_export_revenue_new_exporters5yr.(place{cc})=[];
    share_export_output_new_exporters5yr.(place{cc})=[];
    
    avg_markups_produced_vec.(place{cc})=[];
    avg_markups_consumed_vec.(place{cc})=[];
    dispersion_markups_consumed_vec.(place{cc})=[];
    
    wage_level_vec.(place{cc})=[];
    real_wage_growth_vec.(place{cc})=[];
    real_consumption_vec.(place{cc})=[];
    relative_producer_wage_vec.(place{cc})=[];
    tornqvist_TFP_growth_vec.(place{cc})=[];
end
trade_elasticity_avg_vec = [];
trade_elasticity_US_vec = [];
trade_elasticity_OECD_vec = [];
wage.OECD = 1;